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Frequency  shifts  of  the  199Hg+  5dl06s  2Sm 
(F  =  0 ,Mf=  0)  to  5d96s2  2D5/2  (F  =  2, 

Mf  =  0)  electric-quadrupole  transition  at 
282  nm  due  to  external  fields  are  calcu¬ 
lated,  based  on  a  combination  of  measured 
atomic  parameters  and  ab  initio  calcula¬ 
tions.  This  transition  is  under  investigation 
as  an  optical  frequency  standard.  The 
perturbations  calculated  are  the  quadratic 
Zeeman  shift,  the  scalar  and  tensor 
quadratic  Stark  shifts,  and  the  interaction 
between  an  external  electric  field  gradi¬ 
ent  and  the  atomic  quadrupole  moment.  The 


quadrupole  shift  is  likely  to  be  the  most 
difficult  to  evaluate  in  a  frequency  standard 
and  may  have  a  magnitude  of  about  1  Hz 
for  a  single  ion  in  a  Paul  trap. 
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1.  Introduction 

It  has  long  been  recognized  that  a  frequency  standard 
could  be  based  on  the  282  nm  transition  between  the 
ground  5dl06s  2Sic  level  and  the  metastable  5d96s2  2D5/2 
level  of  Hg+  [1],  The  lifetime  of  the  upper  level  is  86(3) 
ms  [2],  so  the  ratio  of  the  natural  linewidth  An  to  the 
transition  frequency  v0  is  2  X  1 0  ,5.  (Unless  otherwise 
noted,  all  uncertainties  given  in  this  paper  are  standard 
uncertainties,  i.e.,  one  standard  deviation  estimates.) 
Doppler  broadening  can  be  avoided  if  the  transition  is 
excited  with  two  counter-propagating  photons,  as  origi¬ 
nally  proposed  by  Bender  et  al.  [1]  and  subsequently 
demonstrated  by  Bergquist  et  al.  [3].  However,  optical 
Stark  shifts  are  greatly  reduced  if  the  transition  is  driven 
instead  with  a  single  photon  by  the  electric-quadrupole 
interaction.  In  this  case,  Doppler  broadening  can  be 
eliminated  if  the  ion  is  confined  to  dimensions  much 
less  than  the  optical  wavelength,  as  was  first  demon¬ 
strated  by  Bergquist  et  al.  [4]. 


Recently,  the  ( F  =  0,  MF  =  0)  to  ( F  =  2,  M,  =  0)  hy- 
perfine  component  of  the  199Hg+  5dw6s  2Si/2  tO 
5d96s2  2Ds/2  single-photon  transition  has  been  observed 
with  a  linewidth  of  only  6.7  Hz  by  Rafac  et  al.  [5].  A 
laser  servo-locked  to  this  transition  is  an  extremely  sta¬ 
ble  and  reproducible  frequency  reference.  New  develop¬ 
ments  in  optical  frequency  metrology  [6,  7]  may  soon 
make  this  system  practical  as  an  atomic  frequency  stan¬ 
dard  or  clock. 

While  the  (F  =  0,  MF  =  0)  to  {F  =  2,  MF  =  0)  hyper- 
fine  component  has  no  linear  Zeeman  shift,  it  does  have 
a  quadratic  Zeeman  shift  that  must  be  accounted  for.  In 
addition,  there  is  a  second-order  Stark  shift  and  a  shift 
due  to  the  interaction  between  the  electric-field  gradient 
and  the  atomic  electric-quadrupole  moment.  None  of 
these  shifts  has  yet  been  measured  accurately,  so  it  is 
useful  to  have  calculated  values,  even  if  they  are  not  very 
precise.  Also,  it  is  useful  to  know  the  functional  form  of 
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the  perturbation,  even  if  the  magnitude  is  uncertain.  For 
example,  the  quadrupole  shift  can  be  eliminated  by  aver¬ 
aging  the  transition  frequency  over  three  mutually  or¬ 
thogonal  magnetic-field  orientations,  independent  of  the 
orientation  of  the  electric-field  gradient. 

2.  Methods  and  Notation 

The  quadratic  Zeeman  shift  can  be  calculated  if  the 
hyperfine  constants  and  electronic  and  nuclear  g  -factors 
are  known.  Similarly,  the  quadratic  Stark  effect  can  be 
calculated  from  a  knowledge  of  the  electric-dipole  oscil¬ 
lator  strengths.  The  quadrupole  shift  depends  on  the 
atomic  wavefunctions.  Some  of  these  parameters  have 
been  measured,  such  as  the  hyperfine  constants  and 
some  of  the  oscillator  strengths.  There  are  also  pub¬ 
lished  calculations  for  some  of  the  oscillator  strengths. 

Here,  we  estimate,  by  the  use  of  the  Cowan  atomic- 
structure  codes,  values  for  parameters  for  which  there 
are  neither  measured  values  nor  published  calculations. 
The  Cowan  codes  are  based  on  the  Hartree-Fock  ap¬ 
proximation  with  some  relativistic  corrections  [8].  The 
odd-parity  configurations  included  in  the  calculation 
were  5 dwnp  (n  =  6, 7, 8,9),  5dl05f,  5d96s6p ,  5d96s7p, 
5d96s5f,  and  5ds6s26p .  The  even-parity  configurations 
were  5 dwns  (n  =  6,7,8,9,10),  5 dlond  (n  =  6,7,8, 9), 
5d96s 2,  5d96s7s,  5d96s6d,  and  5d96p2.  Recently,  San- 
sonetti  and  Reader  have  made  new  measurements  of  the 
spectrum  of  Hg+  and  classified  many  new  lines  [9]. 
They  also  carried  out  a  least-squares  adjustment  of  the 
energy  parameters  that  enter  the  Cowan-code  calcula¬ 
tions  in  order  to  match  the  observed  energy  levels.  We 
use  these  adjusted  parameters  in  our  Cowan-code  calcu¬ 
lations. 

As  one  test  of  this  method  of  calculation,  we  esti¬ 
mated  the  weakly  allowed  10.7  pm  5dl06p  2Vm  to 
5d96s2  2D3/2  electric-dipole  decay  rate.  This  decay  is  al¬ 
lowed  only  because  of  configuration  mixing,  since  it 
requires  two  electrons  to  change  orbitals.  The  calcula¬ 
tion  shows  the  decay  to  be  due  mostly  to  mixing  between 
the  5dw6p  and  5d96s6p  configurations.  The  calculated 
rate  is  111  s-1;  the  measured  rate  is  52(16)  s-1  [2]. 
Another  test  is  the  electric-quadrupole  decay  rate  of  the 
5d96s2  2D5/2  level  to  the  ground  level.  The  calculated  rate 
is  12.6  s-1,  and  the  measured  rate  is  1 1. 6(0.4)  s-1.  Sim¬ 
ilar  calculations  have  been  carried  out  by  Wilson  [10]. 

Let  H0  be  the  atomic  Hamiltonian,  exclusive  of  the 
hyperfine  and  external  field  effects,  which  are  treated  as 
perturbations.  For  convenience,  we  denote  the  eigen¬ 
states  of  H0  corresponding  to  the  electronic  levels 
5c/ln6.s'  2Sic  and  5d96s 2  2D5/2  having  J:  eigenvalue  M,  by 
IS  1/2  Mj)  and  ID  5/2  Mj ),  respectively. 


The  corresponding  eigenvalues  of  H0  are  denoted 
W(S,  1/2)  and  W(D,  5/2).  An  arbitrary  eigenstate  of  H0 
with  eigenvalue  W(y,J)  and  electronic  angular  mo¬ 
mentum /is  denoted  \yJMj).  Since  199Hg+  has  in  addi¬ 
tion  a  nuclear  angular  momentum  /,  where  I  =  1/2,  the 
complete  state  designation  is  I  y  JFMf),  where  F  is  the 
total  angular  momentum,  and  M,  is  the  eigenvalue  of  Fz . 

3.  Quadratic  Zeeman  Shift 

In  order  to  calculate  the  energy  shifts  due  to  the 
hyperfine  interaction  and  to  an  external  magnetic  field 
B=Bz,  we  define  effective  Hamiltonian  operators 
Hs  and  //rj  that  operate  within  the  subspaces  of  hyper¬ 
fine  sublevels  associated  with  the  electronic  levels 
5dl06s  2Sic  and  5d96s2  2D5/2,  respectively: 

Hs  =  hAsI'J  +  gj(  S)/j.b  J'B  +  g'lp^I-B ,  (1) 

H'd  =  hADI-J  +  gj(  D)/j,b  J'B  +  gipvI'B ,  (2) 

where  As  and  AD  are  the  dipole  hyperfine  constants, 
gj( S)  and  g,/(D)  are  the  electronic  g-factors,  gl  is  the 
nuclear  g  -factor,  h  is  the  Planck  constant,  and  pB  is  the 
Bohr  magneton.  All  of  the  parameters  entering  Hs  and 
III-,  are  known  from  experiments,  although  a  more  accu¬ 
rate  measurement  of  g.,(D)  would  be  useful.  The 
ground-state  hyperfine  constant  As  has  been  measured 
in  a  199Hg+  microwave  frequency  standard  to  be  40 
507.347  996  841  59  (43)  MHz  [11],  The  excited-state 
hyperfine  constant  AD  has  been  measured  recently  by  an 
extension  to  the  work  described  in  Ref.  [5],  in  which  the 
difference  in  the  frequencies  of  the  IS  1/2  0  0)  to  ID  5/2 
2  0)  and  the  IS  1/2  0  0)  to  ID  5/2  3  0)  transition  frequen¬ 
cies  was  determined  to  be  3AD  =  2  958.57(12)  MHz 
1 12],  in  good  agreement  with  an  earlier,  less  precise 
measurement  by  Fabry-Perot  spectroscopy  [13].  The 
ground-state  electronic  g  -factor  g,,(S)  was  measured  in 
198Hg+  by  rf-optical  double  resonance  to  be  2.003  174 
5(74)  [14].  The  excited-state  electronic  g -factor'  gy(D) 
was  measured  in  198Hg+  by  conventional  grating  spec¬ 
troscopy  of  the  398  nm  5dw6p  2P.y2  to  5d96s 2  2D5/2  line 
to  be  1.198  0(7)  [15].  The  difference  in  gj(S)  or  gy(D) 
between  198Hg+  and  199Hg+  is  estimated  to  be  much  less 
than  the  experimental  uncertainties.  The  nuclear  g -fac¬ 
tor  gjis  —5.422  967(9)  X  10^4  [16].  The  measurement 
was  made  with  neutral  ground-state  199Hg  atoms,  so  the 
diamagnetic  shielding  factor  will  be  slightly  different 
from  that  in  the  ion.  However,  this  is  effect  is  negligible, 
since  the  magnitude  of  gl  is  so  small  compared  to  g.,(S) 
or  g./(D). 
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The  determination  of  g.,(D)  could  be  improved  by 
measuring  the  optical-frequency  difference  between 
two  components  of  the  282  nm  line  and  the  frequency  of 
a  ground-state  microwave  transition  at  the  same  mag¬ 
netic  field.  Since  the  uncertainty  in  the  quadratic  Zee- 
man  shift  is  due  mainly  to  the  uncertainty  in  g.,(D),  it  is 
useful  to  see  how  accurately  it  can  be  estimated  theoret¬ 
ically.  The  Lande  g -factor  for  a  2D5/2  state,  including  the 
correction  for  the  anomalous  magnetic  moment  of  the 
electron,  is  1.200  464.  The  Cowan-code  calculation 
shows  that  the  configuration  mixing  does  not  change 
this  value  by  more  than  about  10~6,  i.e.,  1  in  the  last 
place.  There  are  several  relativistic  and  diamagnetic  cor¬ 
rections  that  modify  gy(D),  one  of  which,  called  the 
Breit-Margenau  correction  by  Abragam  and  Van  Vleck 
[17],  is  proportional  to  the  electron  mean  kinetic  energy. 
The  other  corrections  are  more  difficult  to  calculate. 
The  Cowan-code  result  for  the  mean  kinetic  energy  of 
an  electron  in  the  5 d  orbital  of  the  5d96s2  configuration 
is  T  =  19.32  licit,,  where  It,  is  the  Rydberg  constant. 
Using  this  value,  we  obtain  a  theoretical  value  of  gy(D), 
including  the  Breit-Margenau  correction,  of  1.199  85, 
which  disagrees  with  the  the  experimental  value  by 
1.85  X  1CT3,  which  is  2.6  times  the  estimated  experi¬ 
mental  uncertainty  of  Ref.  [15].  If  we  calculate  gj  ( D ) 
for  neutral  gold,  which  is  isoelectronic  to  Hg+,  by  the 
same  method,  we  obtain  a  value  which  differs  from  the 
accurately  measured  experimental  one  [18]  by 
(7  ±  2)  X  1 0  \  Thus,  the  error  in  the  calculated  value 
for  gy(D)  of  199Hg+  might  be  less  than  1  X  1CT4,  but  it  is 
impossible  to  be  certain  of  this,  since  there  are  uncalcu¬ 
lated  terms.  Measurements  of  the  199Hg+  optical  clock 
frequency  at  different  values  of  the  magnetic  field 
should  result  in  a  better  experimental  value  for  gy(D)  in 
the  near  future. 

For  low  magnetic  fields  (B  less  than  1  mT),  it  is 
sufficient  to  calculate  the  energy  levels  to  second  order 
in  B .  To  this  order  in  B ,  the  energies  of  the  hyperfine- 
Zeeman  sublevels  for  the  ground  electronic  level  are 


W( S,  1/2,  0,  0,  B)  =  W( S,  1/2)  - 


fo(S)  ~  g',]2(i&2 


3hAs 


4MS 


W( S,  1/2,  1,  0,  B)  =  W( S,  1/2)  +  ^ 


(3) 


[g.,(S)  -  g/lVgfl2 
4MS 


(4) 


W( S,  1/2,  1,  ±1,  B)  =  W(S,  1/2)  +  ^ 


+  [gAS)  +  g'A^B 
2 

For  the  5d96s2  2D5/2  level  we  have 

W(D,  5/2,  2,  0,  B)  =  Vk(D,  5/2)  - 

[g2(P)  -  glfulB2 
12/7  Ad 

W( D,  5/2,  2,  ±1,  B)  =  VF(D,  5/2) 


7/7  Ad 


[7g.,(D)  -  gl]fiB B  2[g.;(D)  ~  glfulB2 


27  hAu 


W( D,  5/2,  2,  ±2,  B)  =  W( D,  5/2)  - 


7/z  Ar 


[7g,(D)  -  g'A^B  5[g.;(D)  -  g;]2^B2 


108/;  Ad 


W( D,  5/2,  3,  0,  B)  =  W(D,  5/2)  + 


[g/(D)  -  g;f^B2 
12/!  Ad 


W(D,  5/2,  3,  ±1,5)  =  Vk(D,  5/2)  + 


5/!AD 


5/!Ad 


27/i  Ad 


W(D,  5/2,  3,  ±2,  B )  =  Vk(D,  5/2)  + 


5/!Ad 


W(D,  5/2,  3,  ±3,  B)  =  Vk(D,  5/2)  + 


[5g.;(D)  +  g;']/rBB 
2 


5/iAD 


(5) 


(6) 


7/7Ad 


,  (7) 


,  (8) 


(9) 


[5g./(D)  +  g;]/rBff  |  2[gy(D)  -  (1Q) 


[5g./(D)  +  g;]/rBR  5[gy(D)  -  gjfjjjlA 

3  108/7Ad  ’  ^  ; 


(12) 


Here,  W(y,  J,  F,  MF,  B )  denotes  the  energy  of  the  state 
ly  JFMf),  including  the  effects  of  the  hyperfine  interac¬ 
tion  and  the  magnetic  field. 
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At  a  value  of  B  of  0. 1  mT,  the  quadratic  shift  of  the 
IS  1/2  0  0)  to  ID  5/2  2  0)  transition  (optical  clock  transi¬ 
tion)  is  —189.25(28)  Hz,  where  the  uncertainty  stems 
mainly  from  the  uncertainty  in  the  experimental  value 
of  g,/(D).  In  practice,  the  error  may  be  less  than  this  if 
the  magnetic  field  is  determined  from  the  Zeeman  split¬ 
tings  within  the  ID  5/2  F  M,.  )  sublevels.  The  reason  is 
that  an  error  in  g.,(D)  leads  to  an  error  in  the  value  of  B 
inferred  from  the  Zeeman  splittings,  which  partly  com¬ 
pensates  for  the  g,/(D)  error.  If  instead  we  use  the  calcu¬ 
lated  value  of  g./(D),  the  quadratic  shift  for  B  =  0.1  mT 
is  —189.98  Hz,  where  the  uncertainty  is  difficult  to 
estimate. 


4.  Quadratic  Stark  Shift 

The  theory  of  the  quadratic  Stark  shift  in  free  atoms 
has  been  described  in  detail  by  Angel  and  Sandars  [19]. 
The  Stark  Hamiltonian  is 

He=-hE,  (13) 

where  /jl  is  the  electric-dipole  moment  operator. 


(14) 


aScalar(y,  7) 


8  77  €q  v  I(y/ll/xmlly7,)l2 
3(2./  +  1)  yj  W(y' ,  J')  —  W(y,  7)  ’ 


atensor(7,  j )  =  87Te0 


107(27  —  1) 
3(27+3)(7+  1  )(2  7  +  1) 


1/2 


Jl  1  2]  l(y  7II/T1)lly7,)l2 
1 7  7  7'J  W(y',  7')  -  W(y,  7) ' 


(17) 


The  summations  are  over  all  levels  other  than  Iy7). 
Equations  (16)  and  (17)  can  be  rewritten  in  terms  of  the 
oscillator  strengths  fyj,yT- 


,  47T£0e2^2  V 

^scalarv  y  ■>  J )  —  /  i 


.fyJ. 


y'J' 


Tbcnsorl  y ,  7 )  — 


me  ^.[W(y',  7 ')  -  W(y,  7)] 
477e0e2^2 


,(18) 


We 


307(27-  1  )(2  7  +  1) 
(27+3X7+1) 


x2(-ir 


1  1  2  j  f,,yr 

J  J  J'\[W(y',J’)  -  W(y,J)]2 


,  (19) 


where  me  is  the  electron  mass.  The  tensor  polarizability 
is  zero  for  levels  with  7  <  1,  such  as  the  Hg+  5t/'°6.v  2Si/2 

level. 

For  an  atom  with  nonzero  nuclear  spin  /,  the 
quadratic  Stark  shift  of  the  state  \y  JFMf)  is 


and  E  is  the  applied  external  electric  field.  In  Eq.  (14), 
r,  is  the  position  operator  of  the  i  th  electron,  measured 
relative  to  the  nucleus,  and  the  summation  is  over  all 
electrons. 

First  consider  an  atom  with  zero  nuclear  spin,  such  as 
198Hg+.  To  second  order  in  the  electric  field,  the  Stark 
shifts  of  the  set  of  sublevels  \y  JMj)  depend  on  two 
parameters,  ascaiar(y,  7)  and  ctWsoXy,  7),  called  the 
scalar  and  tensor  polarizabilities.  In  principle,  when 
both  magnetic  and  electric  fields  are  present  but  are  not 
parallel,  the  energy  levels  are  obtained  by  simulta¬ 
neously  diagonalizing  the  hyperfine,  Zeeman,  and  Stark 
Hamiltonians.  In  practice,  the  Zeeman  shifts  are  nor¬ 
mally  much  larger  than  the  Stark  shifts,  so  that  HE  does 
not  affect  the  diagonalization.  In  that  case,  the  energy 
shift  of  the  state  \y  JMj)  due  to  //,  is 


AWE{y,  7,  Mf,  E)=-\  ascaiar(y,  7,  F)E2 


4  C^tensorl  y ,  7 ,  F ) 


[3 Mj  -  F(F+  1)] 
F(2F  -  1) 


(3 E2  -  E2).  (20) 


We  make  the  approximation  that  hyperfine  interaction 
does  not  modify  the  electronic  part  of  the  atomic  wave- 
functions  (the  /  7-coupling  approximation  of  Angel  and 
Sandars  [19]).  This  approximation  is  adequate  for  the 
present  purpose,  which  is  to  evaluate  the  Stark  shift  of 
the  199Hg+  optical  clock  transition.  Obtaining  the  differ¬ 
ential  Stark  shift  between  the  hyperfine  levels  of  the 
ground  state,  which  is  significant  for  the  199Hg+  mi¬ 
crowave  frequency  standard  [11],  requires  going  to  a 
higher  order  of  perturbation  theory  [20].  In  the  /  7-cou¬ 
pling  approximation  [19], 


AWE{y,  7,  Mj,  E)=-\  ascata(y,  J)E2 


OscalaJy,  7,  F)  =  0^(7,  7),  (21) 


>Xy,  7) 


[3M,2- 7(7+1)] 
7(27-  1) 


(3  E2-E2).  (15) 


Treating  H ,  by  second-order  perturbation  theory  leads 
to  the  following  expressions  for  the  polarizabilities  [19]: 


«,ensor(y,7,F)  =  (-l)/+',+f 

F(2F  —  1)(2F+  1)(27  +  3)(27+  1)(7+  1)T/2 
(2F  +  3)(F+  1)7(27-  1) 

x{^  J )■ 
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Equations  (18)  and  (19)  were  used  to  evaluate  the 
polarizabilities  for  the  Hg+  5t/  1(,6v  2Si,2  and  5d%s2  2D5,2 
levels.  For  the  calculation  of  Q?scaiar(S,  1/2),  the  oscillator 
strengths  for  all  electric-dipole  transitions  connecting 
the  5dl06s  configuration  to  the  5 dwnp  ( n  =  6,7,8)  and 
5d96s6p  configurations  were  included.  These  were 
taken  from  the  theoretical  work  of  Brage  et  al.  [21],  The 
final  result  is  ascaiar(S,  l/2)/(4TTe0)  =  2.41  X  1CT24  cm3, 
which  compares  very  well  with  the  value  of 
2.22  X  10”24  cm3  obtained  by  Henderson  et  al.  from  a 
combination  of  experimental  and  calculated  oscillator 
strengths  [22],  For  the  calculations  of  ascaiar(D,  5/2)  and 
atensoJD,  5/2),  the  oscillator  strengths  for  electric-dipole 
transitions  to  the  5 dwnp  ( n  =  6,7,8),  5dl05f,  and 
5d96s6p  configurations  were  taken  from  Brage  et  al. 
[21],  The  oscillator  strengths  for  electric-dipole  transi¬ 
tions  to  the  5d96s7p  and  5d&6s26p  configurations  were 
taken  from  the  Cowan-code  calculations.  The  results 
were  ascaiar(D,  5/2)/(4Tre0)  =  3.77  X  1CT24  cm3  and 
atensoJD,  5/2)/(4tt£0)  =  —0.263  X  10-24  cm3.  Evaluating 
Eq.  (22)  for  F  =  2  and  F  =  3  i n  the  5d96s2  2D5/2  level,  we 
obtain  ate„Sor(D,  5/2,  2)  =  fatens0r(D,  5/2)  and  atens0r(D, 
5/2,  3)  =  a,ensor(D,  5/2). 

The  tensor  polarizability  is  much  smaller  than  the 
scalar  polarizabilities  and  in  any  case  does  not  con¬ 
tribute  if  the  external  electric  field  is  isotropic,  as  is  the 
case  for  the  blackbody  radiation  field.  The  net  shift  of 
the  optical  clock  transition  due  to  the  scalar  polarizabil¬ 
ities  is  ^[ttscaiaJS,  1/2)  —  ascaiar(D,  5/2 )]E2.  In  frequency 
units,  the  shift  is  —1.14  X  10-3  E2  Hz,  where  E  is  ex¬ 
pressed  in  V/cm.  The  error  in  the  coefficient  is  difficult 
to  estimate,  particularly  since  it  is  a  difference  of  two 
quantities  of  about  the  same  size.  However,  the  total 
shifts  are  small  for  typical  experimental  conditions.  If 
the  electric  field  is  time-dependent,  as  for  the  blackbody 
field,  the  mean-squared  value  (E2)  is  taken.  At  a  temper¬ 
ature  of  300  K,  the  shift  of  the  optical  clock  transition 
due  to  the  blackbody  electric  field  is  —0.079  Hz.  The 
mean-squared  blackbody  field  is  proportional  to  the 
fourth  power  of  the  temperature.  For  a  single,  laser- 
cooled  ion  in  a  Paul  trap,  the  mean-squared  trapping 
electric  fields  can  be  made  small  enough  that  the  Stark 
shifts  are  not  likely  to  be  observable  [23]. 

5.  Electric  Quadrupole  Shift 

The  atomic  quadrupole  moment  is  due  to  a  departure 
of  the  electronic  charge  distribution  of  an  atom  from 
spherical  symmetry.  Atomic  quadrupole  moments  were 
first  measured  by  the  shift  in  energy  levels  due  to  an 
applied  electric-field  gradient  in  atomic-beam  reso¬ 
nance  experiments  [24,  25]. 


The  interaction  of  the  atomic  quadrupole  moment 
with  external  electric-field  gradients,  for  example  those 
generated  by  the  electrodes  of  an  ion  trap,  is  analogous 
to  the  interaction  of  a  nuclear  quadrupole  moment  with 
the  electric  field  gradients  due  to  the  atomic  electrons. 
Hence,  we  can  adapt  the  treatment  used  for  the  electric- 
quadrupole  hyperfine  interaction  of  an  atom  [26].  The 
Hamiltonian  describing  the  interaction  of  external  elec¬ 
tric-field  gradients  with  the  atomic  quadrupole  moment 
is 

2 

Hq  =  V£(2)-0(21  =  2(- 1  )qVE™0™,  (23) 

<7=2 

where  \7 E<2)  is  a  tensor  describing  the  gradients  of  the 
external  electric  field  at  the  position  of  the  atom,  and 

0O-) 

is  the  electric-quadrupole  operator  for  the  atom. 
Following  Ref.  [26],  we  define  the  components  of 

V£,2)  as 


V42>=-^, 

2  dz 

(24) 

2  V6a^_  +  Y^ 

~  6  dz  ~  6  -  z’ 

(25) 

VE%=-j^d±E±, 

(26) 

where  E±  =  Ex  ±  i Ev  and  d±  =  ^  ±  i^. 

The  operator  components  (9®  are  defined  in  terms  of 
the  electronic  coordinate  operators  as 

&ts2)  =  ~  |  2  (3^/  “ 

j 

(27) 

®±i  =  ~  eJ^'Zzj(xj±  i yj), 

v  z  i 

(28) 

m=-e^2(xj±iyjf. 

(29) 

where  the  sums  are  taken  over  all  the  electrons.  The 
quadrupole  moment  &(y,J)  of  an  atomic  level  1  yj)  is 
defined  by  the  diagonal  matrix  element  in  the  state  with 
maximum  Mj\ 

0(y,j)  =  {yjj\0^\yjj). 

(30) 

This  is  the  definition  used  by  Angel  et  al.  [24], 

In  order  to  simplify  the  form  of  VEa\  we  make  a 
principal-axis  transformation  as  in  Ref.  [27].  That  is,  we 
express  the  electric  potential  in  the  neighborhood  of  the 
atom  as 
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<P(x',  y',  z')  =  A  [(x'2  +  y'2  -  2 z'2)  +  e(x'2  -  y'2)].  (31) 

The  principal-axis  (primed)  frame  (x',  y',  z')  is  the  one 
in  which  <P  has  the  simple  form  of  Eq.  (31),  while  the 
laboratory  (unprimed)  frame  ( x,y,z )  is  the  one  in 
which  the  magnetic  field  is  oriented  along  the  z  axis. 

The  tensor  components  of  VEa>  in  the  principal-axis 
frame  are  obtained  by  taking  derivatives  of  (lHx' ,  y\  z' )'■ 

VE{ry=-2A,  (32) 

=  0,  (33) 

VE?}=^eA.  (34) 


principal-axis  frame.  These  matrix  elements  are  of  the 
form 

(  y  JF/jl  I  &p’  I  y  JF/jl ) 

=  ]>>$.(«)£>$’(«)  '{yJFm'\0fy\yJFm)',  (38) 

m'm 

=  ( y  JF II  ©(2)ll  y  JF  )2  (■ - 1  f~m' 

m'm 

x(  F  ,  2  f)Ow)flW*(W),  (39) 

\—m  q  mj 

=  JF\\@m\\y  JF) 


In  the  principal-axis  frame,  Hq  has  the  simple  form 
Hq  =  —2A&P'  +  sji  eA^OP'  +  (35) 

As  long  as  the  energy  shifts  due  to  Hq  are  small 
relative  to  the  Zeeman  shifts,  which  is  the  usual  case  in 
practice,  Hq  can  be  treated  as  a  perturbation.  In  that 
case,  it  is  necessary  only  to  evaluate  the  matrix  elements 
of  Hq  that  are  diagonal  in  the  basis  of  states  ly  JFMf), 
where  F  is  the  total  atomic  angular  momentum,  includ¬ 
ing  nuclear  spin  I ,  and  MF  is  the  eigenvalue  of  E,  with 
respect  to  the  laboratory  (not  principal-axis)  frame.  Let 
m  denote  the  set  of  Euler  angles  {a,  [3,  y }  that  takes  the 
principal-axis  frame  to  the  laboratory  frame.  To  be  ex¬ 
plicit,  starting  from  the  principal-axis  frame,  we  rotate 
the  coordinate  system  about  the  z  axis  by  a ,  then  about 
the  new  y  axis  by  [3 ,  and  then  about  the  new  z  axis  by 
y  so  that  the  rotated  coordinate  system  coincides  with 
the  laboratory  coordinate  system.  We  can  set  y  =  0, 
since  the  final  rotation  about  the  laboratory  z  axis, 
which  is  parallel  to  B ,  has  no  effect.  The  states  I  y  JFm ) ' 
defined  in  the  principal-axis  frame  and  the  states 
I  y  JF/jl  )  defined  in  the  laboratory  frame  are  related  by 


2  F 
q  m 


i  D^.(ft»)D^_, 


(w), 


(40) 


K  mm'  nri 


x 


(-\f-»-<i(yJF\\0a>\\yJF)  2  (2K+1) 

Af(w), 


E  2  E  E  E  E  E  E  K 
-  in  q  m/  \ /i  —  fi  n  J\m  —m  n 


(41) 

=  ( - 1  j^-'Ty  JE 1 1  @(2)l  I  y  JE )  (  ^  2 

(42) 


where  Eq.  (39)  follows  from  the  Wigner-Eckart  theo¬ 
rem,  and  Eqs.  (40),  (41),  and  (42)  follow  from  Eqs. 
(4.2.7),  (4.3.2),  and  (3.7.8)  of  Ref.  [28],  respectively. 
The  required  rotation  matrix  elements  are,  from  Eq. 
(4.1.25)  of  Ref.  [28]  (with  correction  of  a  typographical 
error), 


Axf(ftJ)  =  1(3  cos2  /3-1), 


(43) 


fl®fti)  =  a/§  sin2  /3(cos  2a  +  i  sin  2a).  (44) 


ly  JFm)' D^(io)lyJFfx), 


(36) 


The  3 -j  symbol  in  Eq.  (42)  is 


where  ( co )  is  a  rotation  matrix  element  defined  in 
the  passive  representation  [28,  29].  The  inverse  relation 
is 


ly  JFfi)  =  2  Dp!,* (o))\y  JFm)' .  (37) 

m 


(  F  2  F 
\-/r  0  i± 

_ 2[3M2-E(E+1)] _ 

[(2E  +  3)(2E  +  2)(2E  +  1)2E(2E  -  1)]1/2 ' 

(45) 


In  order  to  evaluate  the  diagonal  matrix  elements  of 
Hq  in  the  laboratory  frame,  it  is  necessary  to  evaluate 
matrix  elements  of  the  operators  D,'2’’,  defined  in  the 
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The  diagonal  matrix  elements  of  Hq  in  the  laboratory 
frame  are 

(y  JFMF\HQ\y  JFMf) 

— 2[3 Ml  -  F(F+  l)]A(yJF\\0a>\\yJF) 

~  [(2F  +  3)(2F  +  2)(2F  +  1)2F(2F-  1 )] 1/2 

X  [(3  cos2  /S  —  1)  —  e  sin2  f3( cos2  a  —  sin2  a)].  (46) 

It  is  simple  to  show,  by  directly  integrating  the  angu¬ 
lar  factor  in  square  brackets  in  Eq.  (46),  that  the  average 
value  of  the  diagonal  matrix  elements  of  Hq,  taken  over 
all  possible  orientations  of  the  laboratory  frame  with 
respect  to  the  principal-axis  frame,  is  zero.  This  also 
follows  directly  from  the  fact  that  the  quantity  in  square 
brackets  is  a  linear  combination  of  spherical  harmonics. 
It  is  less  obvious  that  the  average,  taken  over  any  three 
mutually  perpendicular  orientations  of  the  laboratory  z 
quantization  axis,  is  also  zero.  This  result  is  proven  in 
Appendix  A.  This  provides  a  method  for  eliminating  the 
quadrupole  shift  from  the  observed  transition  frequency. 
The  magnetic  field  must  be  oriented  in  three  mutually 
perpendicular  directions  with  respect  to  the  trap  elec¬ 
trodes,  which  are  the  source  of  the  external  quadrupole 
field,  but  with  the  same  magnitude  of  the  magnetic  field. 
The  average  of  the  transition  frequencies  taken  under 
these  three  conditions  does  not  contain  the  quadrupole 
shift. 

The  reduced  matrix  element  in  Eq.  (46)  is,  in  the 
//-coupling  approximation, 

(y(IJ)F\\0t2)\\y(IJ)F) 

=  (-ir(2 F+  1){£  ]  '}(/,  l  J^l0(y,J),  (47) 

where  I  is  included  in  the  state  notation  in  order  to 
specify  the  order  of  coupling  of  /  and  /.  For  the  partic¬ 
ular  case  of  the  199Hg+  5d96s2  2DJ/2  level,  the  reduced 
matrix  elements  are 


(D  5/2  2II@(2,IID  5/2  2)  =  2y 

/l4 

Jy  @(D,  5/2), 

(48) 

(D  5/2  3II@(2iIID  5/2  3)  =  2A 

Jy  @(D,  5/2), 

(49) 

Since  the  Cowan-code  calculation  shows  that  there  is 
very  little  configuration  mixing  in  the  199Hg+  5d96s2  2D5/2 
level,  0(1) ,  5/2)  can  be  reduced  to  a  matrix  element 
involving  only  the  5d  orbital: 


0(D,  5/2)  =  |  <5 d2d5l2,  rrij 
=  5/2l3z2  -  r2\5d2d5/2,  nij  =  5/2),  (50) 

=  |  (5 d,  m,  =  2l3z2  -  r2l5/,  m,  =  2),  (51) 

A-tt 

=  eyj^j-(5d,m,  =  2\Y2,0(e,  </>)\5d,  m,  =  2),  (52) 

/  A  rrr 

=  e^—(5d\r2\5d) 

xf  f  F2*2(0,  <^>)F2>o(0,  <^)T2i2(0,^>)sin  9d9d<p,  (53) 
Jo  Jo 

=  5'<5‘'"-'i5‘i)(-2  02X000)-  (54> 

=  -  y  (5d\r\5d).  (55) 

The  apparent  sign  reversal  in  Eq.  (50)  relative  to  Eqs. 
(27)  and  (30)  is  due  to  the  fact  that  the  quadrupole 
moment  is  due  to  a  single  hole  in  the  otherwise  filled  5d 
shell  rather  than  to  a  single  electron .  According  to  the 
Cowan-code  calculation, 

(5/lr2l5</>  =  2.324  al  =  6.509  X  10"17  cm2,  (56) 

where  a0  is  the  Bohr  radius. 

Since  the  quadrupole  shifts  are  zero  in  the  5dl06s  2Si/2 
level,  the  quadrupole  shift  of  the  199Hg+  optical  clock 
transition  is  due  entirely  to  the  shift  of  the  ID  5/2  2  0) 
state,  and  is  given  by 

<D  5/2  2  0I/7qID  5/2  2  0) 

=  |A@(D,  5/2) [(3  cos2  j3  —  1) 

—  e  sin2  (3 (cos2  a  —  sin2  a)],  (57) 

=  -  |Ae(5/l/-2l5/)[(3  cos2/3-  1) 

—  e  sin2  /3( cos2  a  —  sin2  a)],  (58) 

*  -  3.6  X  10~3hA[(3  cos2  /3  —  1) 

—  e  sin2  f3  (cos2  a  —  sin2  a)]Hz,  (59) 

where  A  is  expressed  in  units  of  V/cm2.  Thus,  for  typical 
values  A  ~  103  V/cm2  and  lelA  1,  the  quadrupole  shift  is 
on  the  order  of  1  Hz. 
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6.  Appendix  A.  Angular  Averaging  of  the 
Quadrupole  Shift 

For  the  purpose  of  describing  the  quadrupole  shift, 
the  orientation  of  the  laboratory  (quantization)  axis  with 
respect  to  the  principal-axis  frame  is  defined  by  the 
angles  [3  and  a .  In  the  principal-axis  coordinate  system, 
a  unit  vector  along  the  laboratory  z  axis  is  defined  in 
terms  of  (3  and  a  by 

z  =  (sin  [3  cos  a,  sin  f3  sin  a,  cos  [3 ).  (60) 

We  wish  to  show  that  the  angular  dependence  of  the 
quadrupole  shift  is  such  that  the  diagonal  matrix  ele¬ 
ments  given  by  Eq.  (46)  average  to  zero,  for  £  along  any 
three  mutually  perpendicular  directions. 

An  arbitrary  set  of  three  mutually  perpendicular  unit 
vectors  e\,  e3,  and  e3  can  be  parameterized  by  the  set  of 
angles  #,</>,  and  if/  in  the  following  way: 

e\  =  (sin  9  cos  <f> ,  sin  9  sin  cf>,  cos  9),  (61) 

€2  =  (cos  </>  cos  9  cos  if/  —  sin  </>  sin  if/,  sin  <f>  cos  9  cos  if/ 

+  cos  cf>  sin  if/,  —sin  9  cos  if/),  (62) 

c3  =  (— cos  <f>  cos  9  sin  if/  —  sin  cf)  cos  if/,  —sin  cf>  cos  9  sin  if/ 
+  cos  cf)  cos  if/,  sin  9  sin  if/).  (63) 

It  can  be  verified  by  direct  computation  that  e,-e,  =  S,f. 

The  quadrupole  shift  can  be  evaluated  for  each  of 
these  three  unit  vectors  substituted  for  z  [Eq.  (60)]  and 
the  average  taken.  First  consider  the  average  of  the  quan¬ 
tity  (3  cos2  j8  —  1)  that  appears  in  Eq.  (46):  We  use  the 
fact  that  cos  (3  is  the  third  component  of  £,  so  the 
average  is: 

(3  cos2  (3  —  1)  =  cos 29  +  sin2#  cos2i f/+  sin2#  sin2)//  —  1, 


(64) 

=  cos2#  4-  sin2#  —  1, 

(65) 

=  0, 

(66) 

for  arbitrary  9,  if/,  and  if/.  Similarly,  the  average  of 
the  other  angle-dependent  term  in  Eq.  (46), 
sin2  f3{ cos2  a  —  sin2  a),  is  calculated  by  making  use  of 
the  fact  that  sin  f3  cos  a  is  the  first  component  of  z,  and 
sin  f3  sin  a  is  the  second: 


+  (cos  </>  cos  9  cos  if/  —  sin  cf>  sin  if/)2 

—  (sin  cf>  cos  9  cos  if/  +  cos  cf>  sin  if/)2 

+  (cos  cf/  cos  9  sin  if/  +  sin  cf>  cos  if/)2 

—  (sin  cf/  cos  9  sin  if/  —  cos  cf/  cos  if/)2],  (67) 

=  0,  (68) 

for  arbitrary  9, 1 />,  and  if/.  Hence,  the  matrix  elements  of 
Hq  given  by  Eq.  (46)  average  to  zero  for  any  three 
mutually  perpendicular  orientations  of  the  laboratory 
quantization  axis. 
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